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■ ABSTRACT 

\Q [ This is the second in a series of papers dedicated to unveil the mass structure and 

■^j- , orbital content of a sample of flattened early- type galaxies in the Coma cluster. The 

■"sj" ' ability of our orbit libraries to reconstruct internal stellar motions and the mass com- 

C position of a typical elliptical in the sample is investigated by means of Monte-Carlo 

■ simulations of isotropic rotator models. The simulations allow a determination of the 

optimal amount of regularization needed in the orbit superpositions. It is shown that 

^ c"| under realistic observational conditions and with the appropriate regularization in- 

O ^ ternal velocity moments can be reconstructed to an accuracy of ss 15 per cent; the 

same accuracy can be achieved for the circular velocity and dark matter fraction. In 
contrast, the flattening of the halo remains unconstrained. Regularized orbit super- 
positions are applied to a first galaxy in our sample, NGC 4807, for which stellar 
kinematical observations extend to 3r c ff. The galaxy seems dark matter dominated 
outside r > 2r c g. Logarithmic dark matter potentials are consistent with the data, as 
' well as NFW-profiles (Navarro, Frenk & White 1996), mimicking logarithmic poten- 

k> , tials over the observationally sampled radial range. In both cases, the derived stellar 

j_j ■ mass-to-light ratio T agrees well with independently obtained mass-to-light ratios from 

stellar population analysis. The achieved accuracy is AT « 0.5. Kinematically, NGC 
4807 is characterized by mild radial anisotropy outside r > 0.5r e tr, becoming isotropic 
towards the center. Our orbit models hint at either a distinct stellar component or 
weak triaxiality in the outer parts of the galaxy. 

Key words: stellar dynamics - galaxies: elliptical and lenticular, cD - galaxies: 
kinematics and dynamics — galaxies: structure 
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1 INTRODUCTION 

Interpreting the stellar kinematical data of ellipticals in 
terms of galaxy structural parameters requires knowing the 
gravitational potential as well as the distribution of stellar 
orbits, which - due to projection effects - is not directly 
given by observations. In equilibrium stellar systems, the 
phase-space distribution function (DF) describing the or- 
bital state depends on phase-space coordinates only through 
the integrals of mo tion, admitted by the act ual potential 
(Jeans theorem; e.g. lBinnev fc Tremainel|l987f) . 

In cases where integrals of motion can be expressed (or 

* E-mail: jthomas@mpe.mpg.de 



approximated) in terms of elementary functions, the DF can 
be parameterized explicitly. Several round as well as a cou- 
ple of flattene d ellipticals have been modelled based on this 



approach (e.g.|Dehnen & Gerhardlll994 


:ICarollo et alJ|l995: 


Oian et alJll995l: iDeionehe et alJll996l: 


Gerhard et alJll99c : 


Emsellem. Deionehe & Baconl Il999l: 1 


Matthias & Gerhard 


1999: IStatler. Deionehe & Smcckcr-Hanel 1999: 


Kronawitter et alJ 12000: Saslia et al.l 1 


2000; Gerhard et al. 


200lll. On the other hand. ISchwarzschild'sl dl979l) 



orbit superposition modelling technique provides 
fully general dynamical models for any axisymmet- 
ric or triaxial potential and has been successfully 
applie d to a g rowing number of early-type galaxies 
(e.g. iRix et al] Il997t ICretton fc van den BoschT Il998l : 
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Cretton. Rix fc de Zeeuwl 120001 IGebhardt et~aTI EoO 



Cappellari et alJl2002l: IVerolme et al.ll2002l: IGebhardt et al 



2004 IValluri. Merritt fc Emselleml l2004t iKrainovic et alJ 
2005). 

Even if the number of dynamical models is steadily in- 
creasing, the only comprehensive investigation of elliptical 
galaxy DFs addressing the question of dark matter in these 
systems is still the basis-fu nction based spher i cal m odel- 
ing of 21 round gala xies by iKronawitter et all i2000t) and 
lOerhard et alJ l|200lf) . To extend the results found there and 
in the handful of studies of individual objects quoted above, 
we started a project aimed to probe a sample of flattened 
early-type galaxies i n the Coma clust er, collecting major 
jMehlert et al.l l200(t) and minor axis JWeener et al] 120021) 
kinematical data. The goal is to investigate the dynamical 
structure and dark matter content of these galaxies. 

Most present-day orbit superpositions, conceptionally 
based on Schwarzschild's original implementation, do not 
automatically provide the entire phase-space DF, but only 
orbital occupation numbers or weights, respectively, charac- 
terizing the total amount of light carried by each orbit. In 
principle, changing the orbit sampling strate gy allows one 
to in fer the DF from any orbit superposition iHafner et alJ 
2000), but this approach has not yet been followed fully. To 
take advantage of both, the full generality of orbit super- 
positions and the availa bility of DFs, we extended th e or- 
bit superpositio n code of iRichstone fc Tremaind <ll988T) and 
IGebhardt et al.l i200(J) to reconstruct pha se-space DFs in 
the axisymmetric case (iThomas et alJl2004) . 

Recovering stellar DFs from photometric and kine- 
matic observations - in particular the application of non- 
parametric methods like orbit superpositions - invokes reg- 
ularizatio n in order to pick up smooth phase-space dis tribu- 
tions fe.g. IRichstone fc Tremainelll98gl : lMerritdll993T) . Dif- 
ferent regularization schemes have been applied in the con- 
text of orbit libraries, among them variants o f minimizing 
occupation number gra dients in orbit space Zhac Il996l) 
and maximum entropy iRichstone fc Tremaine Tl38i ) . The 
proper amount of regulariza tion has thereby c ommonly been 
adjus ted to the data (e.g. iRix et ail 11997k IVerolme et alJ 
2002; Richstone et al., in preparation). 

One aim of this paper is to readdress the question of 
how much regularization is needed to recover galaxy in- 
ternal structural properties from observations with spatial 
sampling and noise typical for our sample of Coma galax- 
ies. To this end we study observationally motivated refer- 
ence models under realistic observational conditions and 
optimize regularization with respect to the reconstruction 
of intrinsic input-model properties (see e.g. iGerhard et alJ 
ll998l:ICretton et al.ll999h . By simulating and recovering ref- 
erence galaxies we also determine to which degree the inter- 
nal mass structure and orbital content are constrained by 
observational data in our sample. Furthermore ^-statistics 
are measured to assign confidence intervals to real galaxy 
orbit superpositions. 

To demonstrate the prospects of regularized orbit su- 
perpositions for our project we also present an application 
to one galaxy. We have chosen the faint giant E2 ellipti- 
cal NGC 4807 {Mb = -20.76 for H = 69km/s/Mpc from 
Hyperleda) for the following reasons: (1) The galaxy has a 
prominent boxy photometric feature in the outer parts con- 



straining its inclination. (2) Stellar kinematic data reach out 
to 3 r c ff along the major axis probing its dark halo outside 
the galaxy main-body. (3) NGC 4807 lacks significant minor- 
axis rotation and isophotal twist and thus is an ideal target 
for axisymmetric modelling. (4) The galaxy is only mildly 
flattened and the dynamical models can be compared with 
earlier studies of similar galaxies done mostly in the spheri- 
cal approximation. 

The paper is organized as follows. Sec.|!2]summarizes the 
observations of NGC 4807, Sec. |3] outlines the orbit super- 
position technique. In Sec. [I] Monte-Carlo simulations per- 
formed to derive the optimal regularization are described 
and the achieved accuracy of galaxy reconstructions fol- 
lows in Sec. |^| Orbit models for NGC 4807 are presented 
in Sees. 16181 The paper is closed with a combined summary 
and discussion of the results in Sec.|5] 



2 NGC 4807: MODEL INPUT 

In the following we assume a distance d — 100 Mpc to NGC 
4807 (corresponding to Ho = 69 km/s/Mpc) and take r c g = 
6.7 arcsec as its effective radius dMehlert et aljEoOOt) . 



2.1 Photometric data 

The photometric data are combined from two different 
sources. For the outer parts of NGC 4 807 Kron-Cousins Rc 
band CCD photometry is drawn from lMehlert et "all (2000), 
consisting of profiles for the surface brightness fin , ellipticity 
e an d isophotal shape paramete rs 04 and a& out to « 5.5 r e fr 
fsee lBender fc Mollenhofll987l for a definition of a n ). A see- 
ing of 2 arcsec during the observations causes the profiles 
to be unresolved in their central parts (at d — 100 Mpc one 
arcsec corresponds to 0.485 kpc). To increase the central res- 
olution, the ground-based data are supplemented by corre- 
sponding profiles for fiv, e, 0,4 and ae extracted from archival 
HST V band data (Principal Investigator: John Lucey; Pro- 
posal ID: 5997). 

The two surface brightness profiles /j,r and fiy are 
joined by shifting the HST V band according to the average 
(fiR — /iv) taken over the region 0.75 r c g ^ R ^ 3 r c s , where 
both data sets overlap and seeing effects are negligible. The 
shift {fin — nv) is well defined with a rms of only 0.015 mag. 

For the orbit models ground-based photometry is used 
outside R 3r c ff and HST data inside R < 3r c g. Fig. 
shows the combined {ir, e, 0,4 and ae profiles applied for the 
modelling. 



2.2 Deprojection 

In the implementation of Schwarzschild's modelling tech- 
nique used here, orbit models are not directly fitted to ob- 
served photometry. Instead, the deprojected luminosity den- 
sity is used as a boundary condition for any orbit superpo- 
sition (see Sec. 13. 411 . 

Deprojection of axisymmetric ga laxies is uniqu e only for 
edge-on systems (inclination i = 90° : iRvbickil 1987T) . For any 
inclination i < 90° disk-like konus-densities can be added to 
the luminosity profil e without affecting its pr ojected appear- 
ance for any i' < i iGerhard fc Binnevlll996h . From i = 90° 
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Figure 1. Photometric data for NGC 4807 (dots with error bars, 
from top to bottom: surface brightness fin, ellipticity e, isophotal 
shape parameters 04 and ag) and three deprojections: edge-on 
deprojection (solid); i = 50°-deprojection (dashed); disky i = 
50°- deprojection (dotted). 



to i = 0° the variety of different luminosity models project- 
ing to the same galaxy image generally increases. 

De projections for NG C 4807 are obtained with the pro- 
gram of iMagorrianl Jl99ST) . The code allows one to explore 
the full range of luminosity densities consistent with the pho- 
tometric data by forcing the deprojection towards different 
internal shapes. For example, at any given inclination i the 
goodness of fit of the deprojection can be penalized towards 
any degree of internal boxiness and diskiness, respectively. 
The deprojections are fitted without seeing-correction, since 
our joint photometry includes ground-based data only out- 
side R ^ 3r c ff, where seeing effects are negligible. Outside 
the last measured photometric data point (R > 5.5r e ff), 
the photometry is extrapolated by a de Vaucouleurs profile 
fitted to the inner parts of the galaxy. The isophotes out- 
side R > 5.5 R c a are assumed to be perfect ellipses with a 
constant flattening corresponding to the galaxy's outermost 
measured ellipticity. 

From the average flattening (q) = 0.8 of NGC 4807 we 
expect an inclination angle i > 38°. Lower viewing angles 
would require density distributions intrinsically flatter than 
E7, which are not observed. In Fig. Q three representative 
deprojections are overplotted to the photometric data: the 
(unique) i — 90°-deprojection; a deprojection at i = 50° 
without any shape penalty; a disky i — 50°-deprojection. 
All three luminosity models are equally consistent with the 
data. The differences between the models are illustrated in 
Fig. [5] which displays the isophotal shape parameters and 
the corresponding isophotes of the deprojections as they ap- 
pear viewed from edge-on. At this viewing angle internal 
density distortions show up the strongest. 

As Fig. [5] reveals, both deprojections at i — 50° are 
heavily boxy around R « 2r e fi (en m — 10 per cent for the 



deprojection without shape preference and 04 w —8 per cent 
for the disky deprojection). These distortions are a reflec- 
tion of the drop in 0,4 at R > 3r e ff in the data of NGC 
4807. Projected density distortions progressively strengthen 
in deprojection from i — 90° to i = 0° and assuming i = 50° 
already causes a considerable amplification. Near the center 
- where the observed isophotes are consistent with being 
purely elliptical - the disky deprojection appears smoother 
with on average larger 04, albeit lacking the strong ei4-peak 
occurring in the deprojection without shape preference at 
R = 0.05 7-cff. 

We have constructed orbit models for both, the i = 90° 
and the non-disky i = 50° deprojections, but the heavily 
distorted density distributions at i = 50° lead us to expect 
that we likely view NGC 4807 close to edge-on. 



2.3 Kinematic data 

The kinematic data are described in iMehlert et ail (l2000l) 
and lWeener et all <l2002l) . They consist of two long-slit spec- 
tra for the major and the minor axis, respectively. Pro- 
files of lo wer order Gau s s-Hermite coefficients 70, v, a , Hz 
and Ha dGerhardl Il993l: Ivan der Marel fc Franxl Il993l) pa- 
rameterizing the line-of-sight velocity distribution (LOSVD) 
reach out to 3r c g along the major-axis and out to 0.6r e ff 
along the minor-axis. Our orbit models are not directly 
fitted to the observed Gauss-Hermite parameters, but to 
binned LOSVDs. Therefore the Gauss-Hermite parameters 
are transformed into suitable LOSVDs as follows. First, the 
Gauss-Hermite series according to the observed 70, v, a, H3 
and H4 is evaluated at 1000 values of line-of-sight velocity 
«ios, linearly spaced in the range —4a ^ vi os ^ 4a. At each 
of the sampled projected velocities error bars ALOSVD are 
assigned to the LOSVD from the errors in the corresponding 
Gauss-Hermite parameters by means of Monte-Carlo simu- 
lations. The LOSVD is then binned appropriately for com- 
parison with the orbit library (see Sec. l3.1l below). 

For NGC 4807, the scatter in kinematical data from 
different sides of the galaxy is mostly negligible compared 
with the statistical errors in the data. Therefore, the or- 
bit models are fitted to a symmetrized data set, derived 
by averaging the Gauss-Hermite parameters from each side. 
The largest scatter appears in the outermost major-axis H4- 
measurement. To infer the impact of its uncertainty on the 
reconstructed DF, we separately refitted the mass model 
that most closely matches the symmetrized data set, to the 
kinematics of the two different sides of NGC 4807. Both fits 
lead to similar models in terms of DF and internal kinemat- 
ical structure, suggesting that the asymmetry in this single 
measurement does not largely affect our results. 



3 ORBIT SUPERPOSITION MODELS 

Our method of setting up orbit libraries to construct best 
fitting models is des cribed in detail in Ric hstone et al. (in 
preparation) and in iThomas et al] <l2004l) . Here, we only 
briefly outline some aspects of the method relevant to the 
present paper. 

To recover the mass distribution of a given galaxy by 
orbit libraries, a grid of (parameterized) trial potentials is 
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grid in radius is logarithmic in the outer parts and becomes 
linear in the inner parts; the angular bins cover equal steps 
in sini? (see Richstone et al, in preparation, for details). For 
NGC 4807 the models are calculated in N r x N$ = 200 bins 
with N r = 20 and N$ — 10. Each spatial bin of the model's 
sky-projection is subdivided into JV vo i = 15 bins in projected 
velocity for the purpose of fitting the library to the observed 
kinematics. 



3.2 Luminous and dark mass distributions 

Luminous matter. We assume that the luminous mass of 
NGC 4807 is sufficiently characterized by a constant stellar 
mass-to- light ratio T (see also Sec. [5] and Fig. 1221 for a ver- 
ification of this assumption). The stellar mass density then 
reads T v, where the luminosity density v is taken from the 
deprojections of Sec. 12.21 

Dark matter. Modelling of spiral and dwarf galaxy 
rotation curves indicates shallow inner dark matter 
density distributions, i.e. logarith mic density slopes 
7 = dln(ppM)/d ln(r) « (e.g. Ide Blok et all 1200 ll : 
iGentile et aH l2004). Flat dark matter density cores also fi t 
to the kinemati cs o f round ellipticals jGerhard et alJl200lfl . 
although iRix et alJ II1997I) found steeper profiles consistent 
with observations of NGC 2434. Some optical or radio rota- 
tion curves also allow densit y slopes up to 7 ^ — 1, but they 
do not prefer steep profiles dSwaters et aljl2003T) . 

Shallow dark matter density profiles conflict with 
predictions of (pure dark matter) cosmological TV-body 
simulations. In the ACDM scenario, dark matter dis- 
tributions with central density cusps ar o und —1.5 < 
7 < —1 are found (e.g. iMoore et alJ Il99st iNavarro et alJ 
2004) . Similar steep profiles emerge in warm dark matte r 
JColin. Avila-Reese fc VaienzuelJl200(i iKnebe et aljl2002T) . 
whereas self-inter acting dark matter offers central 7^0 
JPave et al.ll200ll) . 

To probe the whole diversity of theoretically and ob- 
servationally motivated profiles we allow for two different 
dark matter distributions, one with central 7 = 0, represen- 
tative of shallow profiles and another with central 7 = — 1, 
representative of the steeper cases. 



Figure 2. Projected appearance of the three deprojections shown 
in Fig.^when seen edge-on. The top panel shows isophotal shape 
parameters (cf. Fig.0, the bottom panel displays the correspond- 
ing isophotes. Solid: i = 90° deprojection; dashed: i = 50°- 
deprojection; dotted: disky i = 50° -deprojection. 



probed. In each trial potential (Sec. I.' i . 2fi a large set of or- 
bits is calculated (Sec. 13.31 and an orbit superposition is 
constructed that best matches the observational constraints 
(Sec. Th e best-fitting mass parameters with corre- 

sponding errors then follow from a x 2 - ana lysis- 



3.1 Basic grids 

For comparison with observations the meridional plane of 
the orbit model as well as its projection onto the plane of the 
sky are divided into bins in radius r and polar angle i9. The 



Cored profiles. A dark matter distribution that pro- 
vides asymptotically flat circular velocity curves in combi- 
nation with flat (7 = 0) inner density cores is g iven by the 
logarithmic potential jBinnev fc Tremainejll987 



$log(-R,z) = yln 



9 Z 

•I + R 2 + — 
q 2 



(1) 



where R — rcos(#), 2 = rsin(i?) are cartesian coordinates 
in the meridional plane and 1? = 0° corresponds to the 
equator. The density dist ribution generating $log reads 
jBinnev fc TremainejlT987i) 



Plog(R, z) oc v; 



(2q 



2 + l)rl + R 2 + 2(1 - 



1 

2^ 



(r| + R 2 + z 2 /q 2 



) 2 q 2 



(2) 



and is positive everywhere for q G [1/ \/2, 1] . The flatten- 
ing of the density distribution differs from that of the 
potential q. It is generally smaller and varies with radius. 
In the following, we will only consider cored profiles with 
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q — 1.0 (spherical). Together with the assumption of a con- 
stant mass-to-light ratio T models with dark matter distri- 
butions according to equation (J2J can be seen as analogues 
to maximum-disk models of spiral g alaxies and have bee n 
called maximum-stellar-mass models jGerhard et al]ll998lL 

Cuspy profiles. A representative cuspy mass distri- 
bution fitting simulated dark matter halos over a wide 
range of radii is the Navarro- Frenk-White (NFW) profile 
jNavarro. Frenk fc White!ll99rJ) 



PNFw(r,r s ,c) oc 



(r/r s )(l + r/r s ) 2 ' 



(3) 



with a central logarithmic density slope 7 = —1. The pa- 
rameter S c in equation © is related to a concentration pa- 
rameter c via 



200 



3 ln(l + c) - c/(l + c) ' 



(4) 



By the substitution r — » r y/ cos 2 (■#) + sin 2 ($)/q 2 equation 
(O provides halos with a constant flattening q. 

In CDM cosmology the two parameters c and r s turn 
out to be correlated in the sense that higher mass ha- 
los are less concentrated, w ith some scatter due to differ- 
ent mass assembly histories <|Na varro 1 Frenk fc Whitejll996t 



l.ling fc Sutoll2C 

relation reads 



r?oclO (A - logc)/J 



( 20 ofc> 



The corresponding 



(•>) 



Here, we take A = 1 .0 5 and B = 0.15 

{Navarro. Frenk fc White! Il99fj iRix et alJ Il997t) . which 
- for the concentrations 5 < c < 25 considered here - is 
equivalent to within 10 per cent to the relation given in 
IWechsler et~ail (120021) for the now standard ACDM. 

Total gravitating mass and potential. Luminous and 
dark matter components combine to the total mass density 



p = T v + p D M, 



(6) 



with pDM being either pLOG or pnfw. The potential $ fol- 
lows by integrating Poisson's equation. 



3.3 Orbit collection 

Given a large set of orbits is calculated in order to sam- 
ple the phase-space of the potential. Energies E and angu- 
lar momenta L z of the orbits are chosen to connect every 
pair of equatorial radial grid bins by at least one equatorial 
orbit. The surfaces of section (SOS) connected to pairs of 
(E,L Z ) - here the upward orbital crossings with the equa- 
tor - are densely fil led with orbits of all available shapes 
llThomas et al.ll2004h . A typical orbit library contains be- 
tween 6500 and 9500 orbits for L z > 0. Each orbit's retro- 
grade counterpart with L z < is included in the library by 
reversing the azimuthal velocities appropriately. In total a 
typical library then contains between 13000 to 19000 orbits. 



3.4 Orbit superposition 

Any superposition of a library's orbits generates a model 
with a specific internal density distribution and specific pro- 



jected kinematics that can be compared with the observa- 
tions. The relative contribution of each orbit to the superpo- 
sition - the orbital weight Wi - represents the total amount 
of light carried by orbit i. To fit a li brary to a given dataset 
the m aximum entropy technique of iRichstone fc Tremaine! 
( 1988) is applied by maximizing 



S = S — a xlosvd, 
where 



Xlosvd = J]] J]] 
i=i k=i 



mod dat 
^-"-dat 



-jfc 



(7) 



(8) 



quantifies the deviations between the model LOSVDs £ mo d 
and the observed LOSVDs £dat at Njz spatial positions j and 
in the iV vc i bins of projected velocity fc. By 5*, the Boltzmann 
entropy 



S= f fki(f) d 3 rd 3 v =J2 Wiln {^) 



(9) 



of the library's DF / is denote d, Vj is the orbital p hase- 
volume of orbit i, computed as in lThomas et alJ i2004l . and 
a is a regularization parameter (see Sec.^J. 

The decomposition of the library potential $ into two 
components generated by the stellar and the dark matter 
mass distributions, respectively, is meaningful only if the 
final orbit model self-consistently generates the stellar con- 
tribution to the potential. Therefore, the luminosity density 
v is used as a boundary condition for the maximization of 
equation Q. This also guarantees a perfect match of the 
orbit model to the photometric observations. 



3.5 Comparing model with data kinematics 

Since our original data set for NGC 4807 consists of Gauss- 
Hermite parameters up to H4 we will quote the deviations 
between model and data kinematics in terms of 



2 

Xgh 



Nr. 

E 



JJ 3,mod 



A«dat 



J - 1 3,dat 



+ 



mod dat 



+ 



+ 



H 3 

JJ 4,mod 



H 3 

J - 1 4,dat 



A^i.dat 



(10) 



The sum in equation I10H over Nc LOSVDs contains 
■Wdata = 4 x Nc terms and the parameters u mo d, cr mo d, 
#3, mod and //4,mod are obtained from the corresponding 
LOSVDs of equation © by fitting a 4th-order Gauss- 
Hermite series 1 . Equation I10H can only be applied if the 
sampling of the LOSVD is sufficient to get reliable and un- 
biased estimates of the Gauss-Hermite parameters and it 
should be noticed that in the implementation of the or- 
bit superposition method applied here (in contrast to the 



1 The fifth parameter of the fit, the intensity 70, is not included 
in equation I1UI since we scale data as well as model LOSVDs 
to the same surface brightness before comparison. This does not 
automatically imply that the fitted intensities 70 of model and 
data LOSVDs are identic al, but it largely affect s their differences. 
Following the approach o flGerhard et alj | |l998ft we therefore omit 
70 in Xqh- 
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programs building up on the work of iRix et alJ Il997l and 
ICretton et ail Il999t) 

Xgh is n °t explicitly minimized. We 
discuss the effects of using equation (H0> instead of equation 
© to derive confidence regions in Sec. 19.31 



4 REGULARIZATION 

The regularization parameter a in equation J7J controls 
the relative importance of entropy maximization and \ 2 - 
minimization in the model. Basically, increasing a puts more 
weight on the ^-minimization j n the orbit superposition 
and reduces the influence of S in equation 0. For data sets 
with relatively sparse spatial sampling as considered here, 
the freedom in the orbit superpositions allows to fit models 
to the noise in the data when applying a > 1. Such orbit 
models show, however, large density depressions and con- 
tradict the traditional view of relaxed dynamical systems. 
On the other hand, models with lower a have smoother dis- 
tribution functions in the sense that the adopted form for 

5 tends to isotropize the orbital DF, thereby reducing its 
dependency on L z and 73. Fitting orbit models to data sets 
with different spatial coverage and quality will likely change 
the effect of a on the final fit. Likewise, changing the func- 
tional form of S can be used to bias the m odels towards other 
than isotropic DFs frh omas et alJ2004) . The best choice for 
q and S has to be investigated case-by-case, depending on 
the galaxies under study and on the amount and form of 
information that is to be extracted from the observations. 



4.1 Motivation 

Figs.|3]and|3]exemplify the effects of regularization in terms 
of two distribution functions reconstructed from fits of the 
same library for an edge-on isotropic rotator model of NGC 
4807. The Gauss-Hermite profiles in the lower panels of the 
figures are derived from velocity moments obeying higher- 
order Jeans e quations in the self-consiste nt potential of the 
deprojection iMagorrian fc BinnevllT993) . Before calculat- 
ing the Gauss-Hermite parameters the velocity moments are 
slit averaged and seeing convolved to simulate the obser- 
vations of Sec. 12.31 Noise is added to v and a according 
to the fractional errors of the observations and to H3 and 
H4 according to the absolute errors. The lower panels of 
Figs. |3] and |3] respectively, display the rotator kinematics 
along the major and minor axes, together with fits of an 
orbit library containing 2 x 6922 orbits. For Fig. [3] the fits 
were obtained at a — 0.02 and for Fig. 01 at a = 6.73. The 
distribution functions reconstructed from the two fits are 
plotted against orbital energy E in the upper panels of the 
figures. Each dot represents the phase-space density Wi/Vi 
along a single orbit, scaled according to ^ Wj = ^ Vj = 1 
jThomas et ail I2004T) . The more regularized fit in Fig. 
yields XGH/Wdata = 0.3 (iV data = 48), while for Fig. g] the 
goodness of fit is XGH/^data = 0.17. 

The DF at a = 6.73 has density depressions of several 
orders of magnitude for orbits with roughly the same en- 
ergy. Orbit models with such DFs are difficult to interpret as 
the result of relaxation processes like violent relaxation and 
othe rs, occurring in th e dynamical evolution of real galaxies 
(e.g. lLvnden-Belil967l) . They often predict uncommon kine- 
matics along position angles not covered by observational 
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Figure 3. Reconstructed phase-space densities of individual or- 
bits against energy (upper panel) for an isotropic rotator model 
of NGC 4807 (details in the text). For the underlying fit (lower 
panel; dots/squares: major/minor axis isotropic rotator model; 
solid/dashed lines: major/minor axis orbit model) a regulariza- 
tion of a = 0.02 is used. 



data. To illustrate this, the projected kinematics predicted 
along a diagonal axis (position angle i9 = 45°) by the two 
DFs of Figs.|2land2]are plotted in Fig. EI F° r comparison the 
(undisturbed) kinematics of the isotropic rotator model are 
also overplotted. The profiles of the almost non-regularized 
(a = 6.73) fit show large point-to-point variations, which 
cannot easily be reconciled with the scatter in real obser- 
vations. Also, the mean deviation between the unregular- 
ized orbit superposition and the input isotropic rotator (IR) 
model are larger for the non-regularized one than for the reg- 
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Figure 4. As Fig. [3] but for regularization with a = 6.73. Note 
that the lower boundary of phase-space densities in the upper 
panel is partly due to the program setting all orbital weights 
smaller than u>; < w m i n — 10~ 37 equal to Wi = w mln . 



ularized one, e.g. (H| R — H3} — 0.06 for a = 6.73 compared 
with (Hl R - 77 3 ) = 0.02 for a = 0.02. 

The above example demonstrates the importance of reg- 
ularization in the context of recovering internal dynamics of 
isotropic rotators from sparse data sets typical for our sam- 
ple of Coma galaxies. For NGC 4807 we have determined 
the optimal amount of regularization by means of Monte- 
Carlo simulations of several such isotropic rotator models 
(see Sec. 14.21 below). 

The DFs of isotropic rotator models are of the form 
/ = f(E,L z ) and represent only a minority of all possi- 
ble DFs, since they are constant along 73. The choice of 



Figure 5. Projected kinematics of the two DFs of Figs. 151 and 
|2 along a diagonal axis with position angle 1? = 45°. Solid lines: 
a = 0.02 (-+ Fig. EJ; dashed lines: a = 6.73 Fig. gj; dots: 
kinematics of the isotropic rotator model. Note that the orbit 
superpositions were fitted to the model only along the major and 



minor axes. 



such reference models for galaxies like NGC 4807 (fast ro- 
tating, faint giant elliptical s) is observationally motivated 
jKormendv fc Benderlll996t) and further supported by the 
weak ve locity anisotropy fou nd in previous dynamical stud- 
ies (e.g. iGerhard et al.ll200ll see also Sec. |7]for the case of 
NGC 4807), implying only a mild dependence of the DF on 
73. For more luminous ellipticals (shaped primarily by ve- 
locity anisotropy) or lenticulars embedded in roundish halos 
other reference models should be explored. We will turn to 
this in a future publication. 



4.2 Regularization from isotropic rotator models 

The isotropic rotator models constructed to determine the 
optimal regularization with respect to the analysis of NGC 
4807's internal structure are based on the edge-on as well 
as the i — 50° deprojections. For numerical reasons, how- 
ever, the models at i = 50° are forced to have isophotes 
close to pure ellipses by setting 04 = a§ for the deprojec- 
tion. Six models have been probed, three at each of the in- 
clinations i — 50,90°: one self-consistent, one embedded in 
a LOG-halo and a third embedded in a NFW-halo. Kine- 
matic profiles ar e calculated from higher ord er internal ve- 
locity moments (iMaeorrian fc Binnevl H994F) following the 
procedure described in Sec. 14.11 Orbit libraries are fitted 
to the kinematic profiles and the internal velocity moments 
reconstructed from the fits are compared with the original 
input-moments for various values of a. The optimal balance 
between entropy maximization and x 2 - m i mm i za ti°n is re- 
vealed where the reconstructed internal velocity moments 
are closest to the input model. 

Fig-IUshows the results for a self-consistent model based 
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Figure 6. Match A(a) of internal velocity moments and 
XGn/^data °f the Gauss-Hermite parameters as a function of the 
regularization parameter a. Lines represent the mean over 60 fits 
to Monte-Carlo realizations of the isotropic rotator's kinematic 
profiles; shaded regions encompass 68 per cent of the simulations 
closest to the corresponding minimum in A(o) and Xgh/ '-^datai 
respectively. 



on the i = 90°-deprojection. The upper panel shows the rms 
difference A (a) between original and reconstructed internal 
velocity moments (up to second order) as a function of the 
regularization parameter a. The rms is evaluated between 
the innermost and outermost fitted data points, at all po- 
sition angles #, including intermediate ones, not covered by 
data points. The lower panel illustrates the goodness of fit 
XGH/^data- Solid lines correspond to the mean obtained by 
fitting orbit models to 60 Monte-Carlo realizations of the 
simulated data. Shaded areas comprise the 68 per cent frac- 
tion of best reconstructions (upper panel) and 68 per cent 
fraction of best fitting solutions (lower panel), respectively. 

The best results are obtained for a — 0.02, when inter- 
nal kinematics of the fits follow the input moments to an 
accuracy of about 15 per cent in the mean. Lower a yield 
less accurate reconstructions since the orbit superpositions 
do not fit the data enough (see Xgh m the lower panel). For 
larger a on the other hand, the rms A(a) increases, because 
the library starts to fit the noise in the data. 

On average, the orbit models at a = 0.02 fit with 
XGH/^data < 1. A proper normalization for Xgh is however 
hard to obtain, since we do not exactly know the number of 
free parameters in the library. For a — this number is zero, 
because x 2 does not appear in equation Q while for a — > oo 
there are roughly as many free parameters as orbits in the 
library 2 , because then S becomes negligible in equation Q. 
Thus, we do not try to calculate a reduced x 2 - 



2 We have neglected the boundary condition related to v. This is 
expected to reduce the number of free parameters, but is unlikely 
to break its dependency on a. 



Table 1. Summary of simulations aimed at optimizing the reg- 
ularization parameter a. Columns from left to right: inclina- 
tion of reference model; potential of reference model (SC = self- 
consistent, LOG = logarithmic spherical halo, NFW = spherical 
NFW-halo, details in the text); regularization parameter ao that 
yields the best reconstruction of internal velocity moments; rms 
A(«o) between internal velocity moments of reference models and 
reconstructions achieved with a = «o; goodness of fit XQn( a o)i 
number of orbits N or ^i t used for the modeling. 



Table summarizes the total of simulations done to fix 
a. The first column labels the three models investigated at 
each of the two inclinations i = 50° and i = 90° (SC: self- 
consistent models without halo; LOG: a dark halo with r c — 
5kpc and v c = 265km/s according to equation is added 
to the stars; NFW: dark halo profile like in equation with 
r s = 10.0 kpc, c = 27.0; for all models T = 3.0). The halo 
models are constructed to produce a roughly flat circular 
velocity. In the second column of Tabled the regularization 
Qo is quoted attaining the best reconstruction of the internal 
velocity moments, quantified in the third column as A(ao) = 
min(A(a)). The corresponding goodness of fit XGti( a o) is 
given in the fourth column of the table. The number of orbits 
in each of the fitted libraries is found in the last column. 

On average, the six probed models yield (ao) = 0.0179± 
0.0035. The low scatter seems to indicate that ao does not 
depend strongly on the potentials tested. In the implemen- 
tation of the maximum entropy technique used here, a is 
increased iteratively in discrete steps (see Richstone et al., 
in preparation). For the models of NGC 4807 we apply 
a = 0.0199 « 0.02, which is the closest larger neighbor of 
(ao) in these iterations. 

In the remainder of the paper, by quoting 68 (90, 95) 
per cent confidence levels, we always refer to all orbit mod- 
els whose XGH/-W<iata are below the corresponding maximum 
X 2 -level of the 68 (90, 95) per cent best-matching fits in 
the simulations. For these fits parameters like T, i, r c , v c 
. . . are set equal to the true input values, but in the analysis 
of real galaxies they are varied. Strictly speaking then, the 
X 2 -statistics from the simulations are not directly applica- 
ble to real galaxies. For example, to determine the correct 
statistics for the case where T, r c and v c are varied to find 
the best-fitting mass-model, about 10 2 orbital fits for each 
triple of (Y,r c ,t> c ) are necessary to yield the corresponding 
three dimensional x 2 -distribution. Such simulations however 
cannot be performed in a reasonable amount of time with 
present-day computer power. 

Probably, our confidence regions derived as described 
above overestimate the error budget. For example, in 
the simulations we find 68 per cent of all orbits within 
AxGH/-^data = 0.38 from the mean best-fit values. In con- 
trast, applying classical Ax 2 -statistics for a two parame- 
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ter fit, yields about 95 per cent confidence at the same 
Axcm/^data = 0.38. 



5 RECOVERING ISOTROPIC ROTATOR 
MODELS 

The aim of this Section is to quantify in a practical sense 
to which degree a sparse data set like the one described in 
Sec. 12.31 constrains the internal kinematics and mass struc- 
ture of a typical galaxy in our sample. 

5.1 Mass-to-light ratio T and inclination i 

First, we disregard the possible presence of a dark halo and 
try to recover the mass-to-light ratio and inclination of the 
self-consistent edge-on model in TableQ Therefore, Fig. pre- 
sumes the results of fitting libraries with mass-to-light ratios 
in the range 2.0 ^ T sj 4.0 and inclinations i = 50, 70, 90° 
to the kinematics of this isotropic rotator model. The good- 
ness of fit XGH/^data for each pair (T,i) is averaged over 10 
Monte-Carlo realizations of the kinematic profiles and eval- 
uated at two different regularizations. For the upper panel 
q = 0.02 according to the simulations described in Sec. 14.21 
Horizontal dotted lines in the plot correspond to 68, 90 and 
95 per cent confidence levels as derived from the statistics 
of the simulations in Sec. 14.21 

One can read from the plot that the mass-to-light ratio 
is well recovered: the orbital fits reveal T = 3.0 ± 0.5, where 
the input model has Tir = 3.0. The minimum Xgh occurs at 
the true value independent of the assumed inclination. The 
latter is only weakly constrained. For T = 3.0 an inclination 
of i = 70° results in nearly the same XGH/^data as the fit 
with the input irrt = 90° . Even the orbit models with i = 50° 
can only be rejected with less than 90 per cent confidence. 

For the lower panel the fits are evaluated at a = 0.44 
(this choice of a is motivated in Secl^J. Although at this 
less restrictive regularization the confidence intervals shrink 
as compared with a = 0.02 the mass-to-light ratio is now 
less constrained, T = 3.0 ± 1.0. The same holds for the 
inclination: at T = 3.0 all three probed inclinations are now 
within the 68 per cent confidence interval. 

5.2 Internal kinematics 

To give an impression about the ability to recover internal 
motions (assuming the potential is known) Fig. |H| compares 
internal velocity moments reconstructed from libraries fit- 
ted to the self-consistent edge-on isotropic rotator model 
of Sec. 14.11 with the corresponding moments of the input 
model. Solid lines portray the isotropic rotator model and 
points show the average reconstructed moments from fits to 
60 Monte-Carlo simulations of the kinematic profiles; error 
bars indicate 1 a deviations from the mean. For the recon- 
structions a = 0.02 is used. 

The three upper rows demonstrate that for the internal 
major axis the second order moments or, o z and a v are 
accurately reproduced by the fit over the region between the 
pointed vertical lines indicating the radii of the innermost 
and outermost kinematical points included in the fit. The 
fractional errors of or and o v are below 3 per cent and of o z 
are below 6 per cent. For the internal rotation velocity v v the 



. a= 0.0199 


1 

i=90 _ 




i = 70 


\ 

\ 

- v \ 

x \ \ 
\\ \ 


i=50 - 








. a = 0.4410 


- i=90 . 




i = 70 


\ 

\ 

\ 

\ 

V 


i = 50 




i i , i - 


2 


3 4 


stellar mass 


-to— light ratio T 



Figure 7. Xqh P er data point (N^ta, = 48) as function of stellar 
M/L for three different inclinations i. Top panel: with optimal 
smoothing, bottom panel: weak smoothing. The input 2I-model 
has M/L = 3.0 and is edge-on. Pointed horizontal lines in each 
panel represent the Xqh values enclosing 68, 90 and 95 per cent 
of the Monte-Carlo simulations, respectively. 

fractional errors are larger than 10 per cent at some points. 
The match of the rotation velocity can be increased when 
going to larger a, but then the second order moments start 
to wiggle around the input moments and the overall match of 
reconstructed to original moments becomes worse. Fig. [S] is 
representative also for the remaining internal position angles 
of the libraries. It follows that most of the rms in Fig. HJ 
results from a mismatch in v v . 

The last row shows a comparison of the anisotropy pa- 
rameter (3$ — 1 — ct^ / of, , which vanishes for isotropic rotator 
distribution functions / = f(E, L z ). The reconstructions are 
consistent with /3$ = given the scatter caused by the noisy 
data. The small offset of Af3$ = 0.05 is due to a slight over- 
estimation of the radial velocity dispersion and emphasizes 
how sensitive 0-a responds to small inaccuracies in the ve- 
locity dispersions. 

5.3 Mass distribution 

The next step is to recover the structure of the isotropic rota- 
tor model in the second row of Table Q where a logarithmic 
dark halo is present. To this end, we fitted libraries with dif- 
ferent dark halos to pseudo data sets of the rotator model, 
keeping the stellar mass-to-light ratio constant. Fig. |3] com- 
bines the results of the simulations. It shows, from top to 
bottom, cumulative mass-to-light ratio M(r)/L(r), circular 
velocity t> c irc(r) and dark matter fraction MoM{r)/M(r) as 
a function of radius. Vertical dotted lines mark the bound- 
aries of the spatial region covered with kinematic data, thick 
lines display the input model. The shaded areas have been 
constructed by determining at each radius the minimum and 
maximum of M/L, v c [ IC and Mom/M of all libraries within 
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Figure 8. Internal velocity moments along the major axis of li- 
braries fitted to the self-consistent edge-on isotropic rotator model 
of Sec. l4.lj Dots: mean from fits to 60 Monte-Carlo realizations of 
the kinematic data; shaded regions: 1 a deviation from the mean; 
solid lines: velocity moments of the input model; dotted lines: 
boundaries of the radial region included in the fit. 



the 68 per cent confidence range derived from the simula- 
tions of Sec. 14.21 Libraries with LOG-halos as well as NFW- 
halos have been tried (see below) and for each library fits 
to 10 realizations of the isotropic rotator kinematics were 
averaged. 

The figure demonstrates that in the region covered by 
the data the mass structure of the input model is well re- 
produced, with progressively larger scatter towards the outer 
edge of the data. The uncertainty of the mass-to-light ratio 
is AM/L « 1.0 at lr e ft, the circular velocity is accurate 
to 10 per cent at the same central distance. In the outer 
parts orbit superpositions with larger masses than the input 
model are better consistent with the data than smaller mass 
models. Nevertheless, the dark matter fraction is determined 
to A(Mum/M) = 0.2 even at the outermost data point. 

The top panel of Fig. 1101 displays 68, 90 and 95 per 
cent confidence intervals for the two parameters r c and v c of 
LOG-halos (cf. equation^. Each dot marks a pair of (r c , v c ) 
probed by fitting a library to 10 realizations of the pseudo- 
data as described above. The input model's dark halo pa- 
rameters are marked by the asterix, the circle designates the 
best-fitting orbit model. 

As expected, the halo parameters are not well con- 
strained. Although the best-fitting parameter pair is close 
to the input model, the 68 per cent confidence contour com- 
prises a large set of libraries and remains open to the up- 
per right edge of the plot. This follows fro m a degeneracy 
betw een the two parameters r c and v c (e.g. iGerhard et alJ 
1998). Increasing r c and v c appropriately puts more mass 
in the outer parts of the halos while rendering the regions 
covered by kinematics roughly unchanged. The increasing 




Figure 9. Accuracy of the mass reconstruction. Thick lines: the 
isotropic rotator model marked by the asterix in Fig. 1101 shaded 
areas: 68 per cent confidence regions of fits to pseudo data sets of 
the model. For each mass-model the goodness-of-fit is averaged 
over 10 fits to MC-realizations of pseudo-data. 

width of the shaded areas in the upper two panels of Fig. [5] 
is an illustration of this degeneracy. 

We also fitted one parameter NFW-halos according to 
equation JSJ to the same isotropic rotator model in order 
to examine whether these profiles can be excluded by the 
data set at hand. The lower panel of Fig. 1101 displays the 
results in terms of XGii/-^data as a function of concentration 
index c. Horizontal dotted lines correspond to 68, 90 and 
95 per cent confidence levels, the solid line shows fits with 
spherical halos (q = 1.0), while for the dashed line q = 
0.7. Apparently, neither the halo flattening nor the central 
halo density slope are constrained by the data, because both 
spherical as well as flattened NFW-halos exist that provide 
fits equally as good as the LOG-halos. These NFW-halos 
join smoothly to the mass-distribution recovered in Fig. [5] 
and mimic LOG-halos over the limited spatial region covered 
by kinematic data. 

Fig. 1101 suggests that it is not possible to discriminate 
between LOG and NFW halos insofar as they provide similar 
mass distributions over the spatial region sampled by kine- 
matic data. More extended data sets are likely to reduce this 
degeneracy, but will probably not completely remove it. On 
the other hand, as Fig. UUreveals, the mass distribution and 
composition of the input model can be recovered well, inde- 
pendently of the specific parameterization chosen to emulate 
it. 



6 DARK MATTER IN NGC 4807 

Now, we turn to the analysis of the distribution of luminous 
and dark mass in NGC 4807. Fig. [ill shows the minimum 
XGrr/^data (scaled, see below) at each modelled stellar mass- 
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Figure 10. Recovering an isotropic rotator model's dark halo: r c 
and v c of the input model model are marked by the asterix in the 
top panel, the best-fit reconstruction is marked by the circle. Lines 
in the top panel display 68, 90 and 95 per cent confidence contours 
of the fits. The lower panel shows — for the same input model - 
X 2 of NFW-type fits (solid/dashed: q = 1.0/ q = 0.7); horizontal 
lines: 68, 90 and 95 per cent confidence levels. Each library has 
been fitted to fO realizations of the input-model kinematics. 



Figure 11. Deviations between the observed kinematics of NGC 
4807 and the fitted libraries in terms of Xqh as a f unc tion of T. 
At each T the minimum min(xQ H ) is plotted, separately for the 
case of LOG-halos (thick solid line), one-parameter NFW-halos 
(thick dashed) and for the self-consistent case (thin solid). For the 
top row the regularization parameter is a = 0.02, for the bottom 
row a = 0.45. 



to-light ratio T G {1.0,2.0,2.5,3.0,3.5,4.0,5.0}. Horizontal 
lines indicate the 68, 90 and 95 per cent confidence limits 
derived from the simulations of Sec. 14. 21 For the upper panel 
q = 0.02 is used and LOG-potentials with stellar mass-to- 
light ratios 2.5 ^ T ^ 3.0 fit best to the data. Projected 
kinematics of the best-fitting orbit model are plotted in the 
upper panel of Fig. 1121 together with the data. The model 
is edge-on, has stellar T = 3.0 and a spherical dark halo 
of LOG-type with core radius r c = 6.8 kpc and asymptotic 
circular velocity v c — 300km/s. It fits the data very well 
with (unsealed) min(xGH/Wdata) = 0.17, even too well when 
compared with the average (xGm/-Wdata) = 0.71 expected 
from the Monte-Carlo simulations. Since both, the scatter 
in the kinematic data from different sides of the galaxy (cf. 
Sec. 12.31 and the point-to-point variations in the profiles are 
smaller than the statistical error bars, the latter might be 
slightly overestimated. Hence, in the following all XGH/Wdata 
are rescaled such that min(xGH/^data) = (XGH/^data)- The 
good match of model and data nevertheless reconfirms that 
NGC 4807 is consistent with axisymmetry. 

Applying regularization with a — 0.02, self-consistent 
models (mass follows light) are ruled out by more than 95 
per cent confidence (see thin solid line in the upper panel 
of Fig. lilt . The best- fitting self-consistent model fits with 
(rescaled) XGH/-^data = 1-67 and is compared with the data 
in the bottom panel of Fig. 1121 It has T = 3.5, larger than 
the best-fitting halo model, but fails to reach the measured 
rotation in the outer parts of the galaxy. A drop in the outer 
major axis //3-profile indicates that the model maintains 
the measured dispersion a around the observed v partly by 



putting light on retrograde orbits. This causes the LOSVD 
to fall off too sharply at velocities larger than v and to de- 
velop a wing at negative velocities, resulting in the too neg- 
ative H%. Lack of support for dispersion at large rotation ve- 
locities hints at mass missing in the outer parts of this model. 
Higher mass-to-light ratios, however, would strengthen the 
mismatch between the central regions of the model and the 
innermost (minor axis) velocity dispersions. 

To push the self-consistent model to the same level 
of agreement with the data reached in the upper panel of 
Fig. 1121 regularization must be lowered to a = 0.45. Al- 
though the best-fitting self-consistent model then gives a 
satisfactory fit in the sense of a reasonable XGH/^data, it is 
again ruled out by more than 95 per cent confidence when 
compared with the mean (xGH/-^data} using the Monte- 
Carlo simulations. This, because halo models at a = 0.45 
still give significantly better fits to the data. The only effect 
of lower regularization - as shown in the run of XGH/^data 
versus T in the lower panel of Fig. 1111 - is a slight broaden- 
ing of the minimum region in XGH/^data and, consequently, 
an expansion of the allowed range of mass-to-light ratios, 
yielding now 2.0 ^ T ^ 3.5 as compared with 2.5 T ^ 3.0 
at a = 0.02. Concluding, NGC 4807 cannot be convincingly 
modelled with self-consistent orbit models, even at weak reg- 
ularization. 

In Fig. 1131 (from top to bottom) the 68 per 
cent confidence regions of cumulative mass-to-light ratio 
M(r) / Ln{r), circular velocity v c [ IC and dark matter frac- 
tion Mr>M(r)/M(r) are shown (analogue to Fig. |5J. Within 
1 r e g dark matter is negligible and the dynamical mass- 
to- light ratio equals the stellar one, M(r)/La(r) = T = 
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Figure 12. Comparison of NGC 4807 (dots: major axis kinematic 
data; squares: minor axis) with the overall best-fitting model (up- 
per panel) and best- fitting self-consistent model (lower panel). 
Both orbit models are at i = 90°. For each data point two error 
bars are given: the statistical error from the observational data 
(larger, without marks) and the error bar corresponding to the 
scaled X GH (marks). 



3.0. Between lr c a and 3r e ff dark matter comes up from 
A/dm/A/ = 0.21 ± 0.14 to M OM /M = 0.63 ± 0.13 and com- 
bines with the luminous matter to a roughly flat circular 
velocity curve with « c i rc = 280 ± 30 km/s at lr c ff and 
«circ = 318 ± 48 km/s at the last kinematic data point. 
In the same spatial region the total mass-to-light ratio rises 
from M{r)/L R (r) = 3.8 ± 0.8 to M{r)/L R {r) = 8.4 ± 2.4. 
Dark mass equals luminous mass at roughly two effective 



Figure 13. Mass-structure of NGC 4807. Shaded regions are con- 
structed as for Fig. [5] In the middle panel, the circular velocity 
curve of the best-fitting orbit model (solid line) and its decompo- 
sition into the stellar (dotted) and dark matter part (dashed) are 
displayed. 



radii. Beyond the last kinematic data point, the models are 
not well constrained and the profiles start to diverge. 

In accordance with the isotropic rotator simulations the 
halo parameters r c and v c are not well constrained - demon- 
strated by the 68, 90 and 95 per cent confidence contours in 
Fig. 1141 The tilted contour cones are narrower than expected 
from Fig. 1101 Note however, that we have rescaled the Xgh 
and hence effectively reduced the error bars as compared 
with the simulations. Therefore, also the shaded areas in 
Fig. 1131 are smaller than the corresponding in Fig. El 

In the lower panel of Fig. 1141 the goodness of fit 
XGH"/-^data is plotted versus concentration c for the one pa- 
rameter family of NFW-profiles. As in the simulations of 
Sec. 15.31 the halo flattening is not constrained and with 
a best-fit X<3H/-^data ~ 0.77 spherical as well as flattened 
NFW-halos provide almost indistinguishable good fits to 
the data as the best-fitting LOG-halos. In the innermost 
regions of those libraries that match the observations, we 
find p/pDM > 10; the logarithmic density slope of the best- 
fitting spherical NFW model is 7 = —1.8 at the outermost 
major-axis data point and 7 = —2.41 at the outer edge of 
the orbit library. Hence, as already pointed out in Sec. 15.31 
the mass distributions of the data-allowed NFW models - 
over the modelled spatial region - resemble halos of LOG 
type. Fig. ^| illustrates the data-allowed range of dark mat- 
ter densities along the major axis. Note that the lower limit 
on pdm is likely due to the limited range of profile shapes 
probed in our study. Dark matter distributions with nega- 
tive central density slopes 7 < (not tested here) might also 
provide acceptable orbit models. 

Models at i = 50° are ruled out by more than 68, but 
less than 90 per cent confidence. The best-fitting case is 
confronted with the kinematic observations in Fig. 1161 Rep- 
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Figure 14. As Fig. I1UI but for libraries fitted to the kinematics 
of NGC 4807. The best-fitting LOG-model is designated by the 
ring. In the lower panel results for spherical (solid line) as well as 
flattened NFW halos are illustrated (dashed line: q = 0.7). 




Figure 15. Major-axis stellar mass-density according to the edge- 
on deprojection (solid line, T = 3.0) and data allowed dark matter 
densities. Both, LOG and NFW-halos are included in the plot. 
(The shaded area is constructed analogously to Fig. El) 



resentative for all low-inclination models, it fails to repro- 
duce the minor-axis //4-profile, but the mass structure of 
this model joins to the profiles in Fig. 1181 and the inferred 
structural properties of the galaxy are robust against the 
assumed inclination. 



Figure 16. As Fig. 1121 but for the best-fitting model at i = 50° 



7 STELLAR MOTIONS IN NGC 4807 

Visual inspection of the surfaces of section reveals that none 
of the orbit libraries for NGC 4807 contains a significant 
fraction of chaotic orbits. In the best-fitting models, at all 
inclinations, no indication for chaos is detectable at all. 

The internal orbital structure of the best-fitting edge- 
on halo model near the major and minor axes is shown in 
Fig. 1171 Around the equatorial plane, at radii 1 r c s < r < 
3r e ff, the model is characterized by radial anisotropy. En- 
hanced radial velocity dispersion ay causes /3# w f3 v m 0.3, 
whereas a® w a v . Taking the large rotation velocity into ac- 
count, energy in azimuthal motion turns out to be roughly 
equal to the energy in radial motion, = « 2 ot + a 2 ~ a 2 . 
On the other hand, motion perpendicular to the equator 
is suppressed. Approaching the central parts the velocity 
structure changes to isotropy and inside r < 0.5 r e g stays 
isotropic. 

Close to the minor axis, NGC 4807 appears nearly 
isotropic in the outer parts 1 r e ff < r < 3 r e ff with de- 
creasing radial dispersion towards the center, such that 
(3# w Pip < —0.5 inside r < 0.5r e ff. 

68 per cent confidence intervals analogous to Fig. [5] for 
velocity anisotropies /3$ and /3 V as well as internal rota- 
tion velocities near the major and minor axes are shown 
in Fig. 1181 Most tightly constrained is the major-axis rota- 
tion, rising from the central parts outwards until settling 
constant between r e g (v to t ~ 124 ± 6 km/s) and 3r c ff 
(vrot = 107 ± 4 km/s). Velocity anisotropy is determined 
to at least A/3 ^ 0.2 and the trends in the anisotropy struc- 
ture of the best-fitting model in Fig. 1171 are clearly seen in 
all allowed orbit models. 

Fig. 1191 illustrates the internal close-to-major axis kine- 
matics of the best-fitting orbit model at i = 50° . In contrast 
to the edge-on models j3$ ~ over almost the whole region 
constrained by observations, with a slight raise towards the 
center. Isotropy in the meridional plane (oy w <r#) is ac- 
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Figure 17. Internal kinematics of the best-fitting halo model Figure 18. Confidence regions for meridional velocity anisotropy 

near the major axis (top panel; position angle $ = 2.9°) and near j3$, azimuthal velocity anisotropy j3 v and internal rotational ve- 

the minor axis (bottom panel; i? = 77.1°), respectively. From top locity v Iot (upper panel: near major axis, position angle $ = 2.9°; 

to bottom: velocity dispersions a r (solid), a$ (short-dashed) and lower panel: near minor axis, position angle $ = 77.1°) of NGC 

a v (long-dashed); velocity anisotropy /3#, f} v ; internal rotational 4807. Shaded areas are constructed as for Fig. 151 
velocity v vot = v v . Vertical, dotted lines indicate the innermost 
and outermost radius of kinematic data. Note, that the minor 

axis data reach only out to 0.6r cff . 8 PHASE-SPACE STRUCTURE OF NGC 4807 

In Fig. HIl phase -space densities Wi/Vi of individual orbits 
are plotted against orbital energy E for the best-fitting or- 
bit model of NGC 4807. The top panel only includes or- 
bits whose maximally reached latitude, i? ma x, is lower than 
"#max < 45° . These orbits are confined to a cone with open- 
ing angle i? < 45° around the equatorial plane. The remain- 
companied by suppressed azimuthal dispersion a v . But in ing orbits are plotted separately in the bottom panel. As in 
contrast to the edge-on model of Fig. 1171 the rotation of the Figs. |3] and 2] the phase-space densities are scaled accord- 
actual model causes (v^,)/a^ ~ 1.7, so that this model is ing to Wj = y~) Vj = 1, where Wi and Vi are the orbital 
dominated by azimuthal motion near the equatorial plane. weights and orbital phase volumes, respectively. Whereas 
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Figure 19. Internal close-to-major axis kinematics (position an- 
gle # = 2.9°) of the best-fitting orbit model at i = 50°. Line 
definitions as in Fig. 1171 





>" 



go -10 



-20 



- 



go -10 



-20 1= 




0.00 < 8 < 45.00 



45.00 < 9 < 90.00 



2000 



4000 



-E [km 2 /s 2 ] 



Figure 20. Reconstructed distribution function of the best- 
fitting orbit model at a = 0.02. Each dot represents the phase- 
space density along a single orbit. Upper row: orbits confined to a 
konus with opening-angle i? m ax < 45° around the equator; lower 
row: remaining orbits. 



are plotted against maximum orbital elongation i9 max in 
the meridional plane. The three rows show orbits located 
in three different spatial regions separately. The top row in- 
cludes only orbits whose invariant curves in the surface of 
section are confined to 0.5r c ff < rsos < lr c ff 3 . The spa- 
tial regions for the other two rows are 1 r e g < rsos < 2 r e g 
and 2r c ff < rsos < 4r e ff, respectively. Crosses denote orbits 
with positive L z > while open circles denote orbits with 
negative L z < 0. From the center outwards, orbits with neg- 
ative angular momentum are progressively depopulated as 
compared with their counterparts at positive L z . The differ- 
ence between prograde and retrograde orbits is most notice- 
able near the equatorial plane and disappears completely for 
bulge orbits reaching up to higher latitudes. In the bottom 
panel of Fig. 1211 orbital phase densities are plotted against 
L z /L Z: d IC , L ZjC irc being the angular momentum of a circu- 
lar orbit with the actual energy. Phase-space densities are 
extracted in the same spatial regions as for the upper panel 
of the Figure. It is now apparent that the keel in the DF of 
Fig. I20l is caused by a depopulation of (retrograde) circular 
orbits (at \L Z \ w i 2 , c irc). 

In fact, depopulation of retrograde near-circular orbits 
leading to the just described keel appears in all data-allowed 
DFs, with a slight tendency to be less pronounced in flat- 
tened halos. Even the DFs obtained assuming NGC 4807 
is inclined by i = 50° are featured by this structure, al- 
beit there, the keel is not exclusively made of retrograde 
orbits: also some of the prograde orbits are suppressed. In 
the low inclination DFs a third prominent substructure in 
phase-space emerges in form of a significant fraction of high- 
energy orbits (around E m — 2500 km 2 /s 2 ) being suppressed 
by about 6 orders of magnitude in density, independent from 
L z . This structure partly overlays the keel and is caused by 
the distorted luminosity density at i = 50° since it shows 
up even in the maximum entropy models (a ~ 0) which are 
only forced to reproduce the density profile, but not fitted 
to the kinematics. 

Orbits in the keel together with their prograde counter- 
parts (approximately all orbits i with wf /w~ > 99, where 
W i = w(Ei,±L x ,i, 1 < i < N pro and iVp ro prograde or- 
bits in the model) make only 4.3 per cent of the total model. 
Without these orbits, a slowly rotating bulge (v ss 70km/s 
at r > r c ff) with a slightly peaked velocity distribution 
{H4 > 0.04 at r > r c g) appears (see also Sec. 19. 41 . On the 
other hand, the extracted orbits rotate with about 200km/s 
and produce a narrow, low-dispersion LOSVD. They develop 
a dumbbell-like structure extending to 2 r e ff above the equa- 
torial plane. The fraction of extracted orbits additionally 
obeying i9 max < 10° shrinks to 0.7 per cent of the total 
model. 



the orbits maintaining the bulge of the galaxy follow a DF 
similar to the one shown in Fig. [!|] some of the orbits con- 
fined around the equatorial plane are strongly depopulated. 
Around E w 3000 km 2 /s 2 phase -space densities are up to 20 
orders of magnitude below the main stream DF. 

In the upper panel of Fig. 1211 orbital phase densities 



3 Here we use rsos as a short cut for the radii of or bital equatorial 
crossin gs, defining the surface of section. Fig. 1 in lThomas et all 
( 2004 ) gives an example of a typical surface of section in the actual 
context. The equatorial crossings locate orbits in the meridional 
plane fairly well. 
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Figure 21. The same model as in Fig. E3 Top panel: ph use 
densities versus maximum orbital elongation $ max in the merid- 
ional plane. Bottom panel: phase densities versus "circularity" 
(L z scaled by the angular momentum L z c ; rc of a circular or- 
bit with same energy). Crosses (open circles) denote orbits with 
positive (negative) L z . Orbits are extracted from different spa- 
tial regions of the model as designated in the panels (see text for 
details). For each spatial region the number of prograde orbits 
(crosses) equals the number of retrograde orbits (open circles). 



9 SUMMARY AND DISCUSSION 
9.1 Regularized orbit models 

We have investigated how closely the internal mass dis- 
tribution and kinematic structure of a galaxy can be re- 
covered from a sparse kinematic data set typical for our 
project aimed at investigating a sample of flattened early- 



type galaxies in the Coma cluster. The degree to which or- 
bit models follow internal galaxy properties depends on the 
amount of regularization ap plied in the fits. In th e max - 
imum entropy technique of iRichstone fe Tremaind <ll988l) 
used here, a regularization parameter a controls the relative 
importance of entropy maximization (regularization) on the 
one side, and fit to the data on the other. To find out which 
choice of a is optimal with respect to our goal of determining 
galaxy structural parameters, we simulated observationally 
motivated isotropic rotator models under realistic observa- 
tional conditions. The models are based on a prototypical 
elliptical in our sample, NGC 4807. By varying a in the fits 
of appropriate orbit libraries to pseudo data of the reference 
models, the match between internal velocity moments of in- 
put model and orbital reconstruction can be evaluated as a 
function of regularization. 

Our simulations indicate that the mass-structure of an 
elliptical can be recovered to about 15 per cent accuracy in 
terms of mass-to-light ratio and circular velocity curve, if 
a regularization parameter a = 0.02 is applied. For the dy- 
namical models tested here, the optimal choice of a turns out 
to be roughly independent from the galaxy's gravitational 
potential. The same accuracy as for the mass structure is 
achieved in the reconstruction of internal velocity moments. 
On the other hand, halo flattening and galaxy inclination 
are only weakly constrained by our data. 

Regularization biases orbit models towards some given 
idealized galaxy model, assumed to represent the object un- 
der study reasonably well. For early-type galaxies, regular- 
ization has been commonly implemented to isotropize the 
final fit. How much of regularization is optimal in the re- 
construction of a given galaxy depends on the specific data 
set (spatial coverage, quality) on the one hand and the de- 
gree to which the regularization bias approximates the given 
galaxy on the other. Isotropizing the DF might, for exam- 
ple, be a good recipe for relaxed early-type galaxies but less 
favorable for lenticulars with a significant, dynamically cold, 
component. In that case either the value of a or, preferen- 
tially, the functional form of S has to be reconsidered. In 
any case, and for any regularization technique, the optimal 
balance between fit to data and smoothing of the DF can 
be examined case-by-case from simulations similar to those 
described here. 



ICretton et ail il99Sl> and IVerolme fc de Zeeuwl (12002 

have determined optimal regularization for an implemen- 
tation of Schwarzschild's method with two-integral com- 
ponents by reconstructing similar DFs as used here. The 
resulting regularization has proven to yield plausible re- 
sults in subseq u ent ap p lications (e.g.lCapp ellari et al]|2002j: 
' Verolme et alT l2002t ICopin. Cretton fc Emselleml 120041: 
Cretton fc Emselleml [20041 : iKrainovic et al.ll2005l) . It is dif- 
ficult to quantitatively compare the amount of regulariza- 
tion applied in these works with the results of our simu- 
lations since the regularization techniques differ. However, 
in each case the achieved accuracy in the recovery of test 
model parameters and the derived dynamical structure of 
real galaxies indicate that the applied amount of regulariza- 
tion is comparable. 
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9.2 Luminous and dark matter in NGC 4807 

We then applied our code with the simulation-derived 
regularization to the elliptical NGC 4807. The dynamical 
models require substantial dark matter in the outer parts 
of the galaxy. Evidence for dark matter in the form of flat 
circular velocity curves from (integrated) stellar kin e matic s 
of ellipticals has also been found bv [fy x et alJ il997t) . 
iGerhard et alJ lll99ct). lEmsellem. Deionghe fc B acon 
|l999t) . ICretton. Rix fc, de Zeeuwl feOOOft.lKronawitter et alJ 
fcoooh . ISaelia et all fcOOOT ) and IGerhard et alJ i200lF T^ 
Dynamical models of five galaxies with planetary nebulae 
kinematics complementing stellar kine matics out to ~ 5 r e s 
lead t he authors to different appraisals: iRomanowskv et alJ 
(2003) argue for non-dark matter models consistent with 
three galaxies observed by the plane t ary ne bulae spectro- 
graph, while IPeng. Ford fc Freem an (2004) require dark 
matter in NGC 5128 and the models of NGC 1399 also 
point at a dark halo flSaglia et al.ll20o7iF) . 

S pherical models for the 21 round ellipticals in the sam- 
ple of lKronawitter et al.l (|200ff l reveal dark matter fractions 
of 10 — 40 per cent at 1 r e ff and dark mass equ als luminous 
mass at roughly 2 - 4r eff ijOerhard et al.ll200^ . both com- 
parable to our results for NGC 4807. The halo core density 
p c of our best-fitting model is about 30 per cent below the 
predictions of their p c — Lb relation (taking Mb = —20.76), 
although still consistent with it. Among the allowed LOG- 
halos however, the core densities vary by a factor of ten and 
taking into account NFW-fits, even halos 90 times denser 
then the best-fitting model are consistent with the data (cf. 
Fig. I15H . From the concentration of our best-fitting (spheri- 
cal) NFW-halo nevertheless a relatively low formation red- 
shift Zf « 2.5 follows. Still, we need a larger number of flat- 
tened ellipticals modelled in sufficient generality to recover 
the detailed properties and physical origin of their mass dis- 
tributions. 

To crosscheck our mass decomposition for NGC 
4807 we compared stellar mass-to-light ratios of the or- 
bit superpositions with stellar mass-to-light ratios de- 
termined completely independent. In Fig. 1221 mass-to- 
light ratios from stellar population models dMarastonl 
" 19981 : iThomas. Maraston fc Bended Exil iMarastonl 1200- 



of NGC 4807 ' s maj or-axis spectrum dMehlert et alJ 1200 
iMehlert et al.l 120031) are compared to the best-fitting or- 
bit models (dashed line: stellar mass-to-light ratio in or- 
bit model; hatched region: 68 per cent confidence region of 
surface mass SM over surface brightness SB in orbit mod- 
els; different symbols refer to different initial-mass-functions 
(IMFs) underlying the populations). Systematic uncertain- 
ties stemming from the unknown IMF are roughly compa- 
rable to the statistical errors, indicated for the Kroupa-IMF 
by the pointed area. Only libraries with T = 3.0 are taken 
into account in the figure. 

As the figure demonstrates, the stellar mass-to-light ra- 
tios determined dynamically agree with stellar populations 
following a Kroupa IMF to AT m 0.5. This (1) confirms 
our mass decomposition and (2) justifies a posteriori the 
assumption of a constant stellar mass-to-light ratio. The 
total surface mass of the dynamical models, as indicated 
by the shaded region in Fig. 1221 exceeds the stellar con- 
tribution by far and underlines the evidence for dark mat- 
ter in NGC 4807. Since our best-fitting orbit models have 
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Figure 22. Projected mass-to-light ratio along the major axis of 
NGC 4807. Hatched area: 68 per cent confidence region of surface 
mass (SM) over surface brightness (SB) for all orbit models with 
T = 3.0 (dashed line). For comparison the projected stellar mass- 
to-light ratios from population analysis (see text) are shown for 
two different IMFs: Salpeter-IMF (triangles), Kroupa-IMF (dots). 
The pointed area encompasses the statistically allowed T based on 
the Kroupa-IMF. For the sake of clarity, the similar, but shifted, 
error region for the Salpeter-based T is omitted. 



LOG-halos (or NFW-halos imitating the mass distribution 
of LOG-halos over the kinematically samp led spatial re- 
gion) they have maximum stellar mass (e.g. IGerhard et alJ 
1998). This is in accordance with most previous studies 
that did not require steep central dark matter p rofiles 
jG^rhard^^lJ_^9^jErrrseUerm ^99j; 
Cretton. Rix fc de Zeeuwl l200d: iKronawitter et al.l hood : 



Saglia et alJ I200C and IGerhard et all 1200 ll) . On the other 
hand. iRixetHd I ll997li found such profiles consistent with 
NGC 2434. 



9.3 Comparing the kinematics: xgh versus Xlosvd 

The confidence intervals for the structural properties of 
NGC 4807 have been derived from Xgh as defined in equa- 
tion 1101 . However, our models do not explicitly minimize 
Xgh but instead Xlosvd as given by equation JSJ . In order 
to investigate how the properties of the best-fitting orbit 
model and the corresponding confidence limits of the halo 
parameters depend on the choice of x 2 we have analyzed 
the Monte-Carlo simulations as well as the models for NGC 
4807 also in terms of 



Xlosvd = XLOSvo/^data, 



(11) 



where iV<j a ta is the number of velocity bins for which C dat > 
(cf. equation |SJ. 

Fig. 1231 shows the best- fitting halo parameters and cor- 
responding confidence intervals for the simulated isotropic 
rotator model of Sec. 15.31 The figure is as Fig. 1101 besides 
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Figure 23. As Fig. I1UI but the confidence regions are derived 
from Xlqsvd ( see Sec. 13.51 . 



Figure 24. As Fig. 1141 but the confidence regions are derived 
from Xlosvd ( see Sec. 13.51 . 



that all confidence regions and model-with-data comparisons 
are computed in terms of Xlosvd- Comparing Figs. I10| and 
1231 it turns out that both x 2 -definitions yield identical best- 
fitting models (designated by the rings in the upper panels, 
and given by the minima of the two curves in the lower pan- 
els). Concerning the confidence regions, however, the com- 
putations based on Xlosvd l ea <A to smaller confidence limits 
for both, the logarithmic and the NFW-halos. 

For NGC 4807, the analysis based on Xlosvd leads to 
the results shown in Fig. 1241 Again, besides the different 
X 2 -definition applied, it is as Fig.[H] The best-fitting NFW- 
halos in Fig. l24l are the same as in Fig. 1141 The shapes of the 
Xlosvd curves for q — 1.0 and q = 0.7 in the bottom panel 
of Fig. 1241 prefer halos slightly more concentrated than in 
the XGH _case 5 but, as for the isotropic rotator simulations, 
the differences in the results obtained from Xgh an< A Xlosvd 
are small. 

Regarding the best-fitting logarithmic-halo models 
(marked by the rings in the upper panels) both \ 2 ~ 
calculations give consistent, but, in contrast to the Monte- 
Carlo simulations, not identical results. As already indi- 
cated in the NFW-fits, the logarithmic halo that best fits 
in the sense of Xlosvd ls more concentrated: the dark mat- 
ter fraction inside r e s is Mum/M — 0.35, compared with 
Mum/M = 0.22 for the best-fitting model in the upper 
panel of Fig. 1141 A more striking difference between the up- 
per panels of Fig. I14l and l24l is the closure of the 68 per cent 
XLOSVD" conn c!ence contour to the upper right. The contour 
differences maybe related to the ignorance of 70 in equation 
HOII , since the mismatch between model and data intensities 
progressively increases from the lower left to the upper right 
in the upper panel of Fig. 1141 

Nevertheless, the comparison of the results obtained 
from Xgh an d Xlosvd reveals that both methods - within 
the errors - give the same results. The confidence limits 



quoted in this paper, based on equation 111)1 . are the more 
conservative choice, but maybe slightly too pessimistic. 



9.4 The outer parts of NGC 4807 

Based on a recent refi nement jThomas et all 20041) o f the 
Schwarzschild code o f iRichston^S^Tremame il988tl and 
ICebhardt et"all j200Cft . we recovered a depopulation of ret- 
rograde, near-circular orbits in the phase-space DF of NGC 
4807, giving rise to a keel when plotting orbital phase-space 
densities against orbital energy. The prograde counterparts 
of the depopulated orbits form a dumbbell-like structure ex- 
tending about « 2 r c ff above the equatorial plane. 

To investigate what might cause the depopulation of 
retrograde orbits in the outer parts of NGC 4807, we pro- 
jected a sequence of distribution functions in which all orbits 
involved in the keel (w~[ /w~[ > 99, see Sec. |SJ are repopu- 
lated according to 



i.ow = i x { w i + u>i ) 



(12) 



Since w Ef new + w~ ncvt = wj + u>~ this does not alter the fit 
to the luminosity profile but just levels the relative fraction 
£ of light on the prograde and retrograde specimen of all 
keel-involved orbit pairs. Roughly speaking, equation 1121 
transforms the keel into a narrow cloud of points parallel to 
the main-stream DF with an offset increasing with |£ — 0.5] . 
Figure 1251 shows the resulting projected (major-axis) kine- 
matics for £ = 0.7, 0.94, 1.0. As expected, shifting light from 
prograde to retrograde orbits reduces the amount of rota- 
tion v and lowers H3 in the outer parts, while at the same 
time leading to a larger velocity dispersion a and enhanced 
H4. As shown by the goodness of fit in the lower panel of 
Fig. 1251 the best fit to the data is achieved for £ = 1.0 (all 
retrograde keel-orbits completely depopulated). Comparing 
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Figure 25. Upper panel: projected major-axis kinematics of 
three distribution functions calculated in the gravitational po- 
tential of the best-fitting mass model: (1) all retrograde orbits in 
the keel depopulated completely (solid); relative fraction of light 
on prograde keel-orbit counterparts (2) £ = 0.94 (short-dashed) 
and (3) £ = 0.7 (dotted). Lower panel: goodness of fit (including 
the minor-axis LOSVDs, which are not affected by the redistribu- 
tion of keel-involved orbits) as a function of the relative fraction § 
of light on prograde keel-orbit correspondents; horizontal dotted 
lines indicate 68, 90 and 95 per cent confidence levels. 

the corresponding solid lines in Fig. I25l with the upper panel 
of Fig. I12l reveals that £ = 1.0 provides essentially the same 
fit as the best-fitting orbit model. Consequently, the retro- 
grade keel-orbits can be regarded as completely depopulated 
in our models of NGC 4807. Reducing the relative fraction of 
light on the prograde counterparts of keel-orbits to £ = 0.94 
is only consistent with the observations at the 68 per cent 
confidence level. A further equalization to £ = 0.7 is already 
incompatible with the observed velocity profiles. 

The outer major-axis LOSVDs of NGC 4807 might in- 
dicate weak tria xiality, for example in form of a weak, nearly 
end-on bar (e.g. lBureau fc Athanassoulal20o'3) . A weak bar- 
like structure would also fit to the boxy appearance of 
the galaxy's outer isophotes. Assuming that NGC 4807 is 
slightly triaxial, the keel in the DF might be an artifact of 
the assumption of axisymmetry. Note, however, that NGC 
4807 is consistent with being axisymmetric, albeit then the 
depopulation of single orbit families, especially around the 
equatorial plane, is hard to understand in the course of dy- 
namical processes. 

Depopulation of retrograde orbits in the outer parts of 
the galaxy is accompanied by a change in stellar ages from 
r m 5 Gyr inside r < r e s to t « 10 Gyr at larger radii 
jMehlert et alf2 003l. Towards 3 r e g stellar ages become un- 
certain and no clear trend is visible. It seems also possible, 
that the keel in the DF is an artifact related to a distinct 
stellar component and an associated (slight) change in T. It 
should be noted that because the orbit models are fitted to 
deprojections and stellar kinematics, the derived DFs char- 



acterize the amount of light per volume element in phase- 
space, not the mass-density in phase-space. Accordingly, de- 
pressions in the DF might (at least partly) reflect enhanced 
T. 

A distinct stellar component in the outer parts must not 
necessarily take the form of a bar, but could also be a faint 
axisymmetric stellar disk. Indeed, the photometry as well as 
the orbit models suggest that NGC 4807 is close to edge-on, 
but a dynamically cold, outer stellar disk made of only 0.7 
per cent prograde orbits at latitudes ■& < 10° (see Sec. |HJ is 
consistent with the non-disky isophotes (maxfl4 < 1). A disk 
carrying 2 per cent of the total luminosity would show up 
with max cu « 4 when seen exactly edg e-on and max 0,4 < 1 
for such a disk would require i < 75° jRix fc Whitelll990l) . 
From our simulations it seems impossible to distinguish kine- 
matically between i — 70° and i = 90° and, consequently, 
there is some freedom for an outer disk in NGC 4807. 

In any case further investigations are necessary to ver- 
ify the significance of the keel-structure in the DF of NGC 
4807. (1) Kinematic measurements along intermediate posi- 
tion angles in the galaxy would likely put more constraints 
on the DF. (2) Detailed comparison of orbit superpositions 
with realistic galaxy models are required to estimate how 
relics of triaxiality, faint disks and multi-component (multi- 
T) structures show up in our axisymmetric models. These 
questions will be addressed in a forthcoming publication. 



9.5 Internal stellar kinematics of NGC 4807 

Along the major-axis, NGC 4807 is mildly radially 
anisotropic. Ra dial anisotropy has b e en found in a number 
of e l lipticals l|Gerhard et alJ Il998t l^i^iias^^GeriLarfl 
| l999UCretton. Rix fc de ZeeuJbooal : ICebhardt et alJl200rT 
ISaglia et all l200Cfl . The amount of anisotropy ((3 m 0.3) 
in the outer parts as well as the isotropy of the central 
region (inside 0.3 r e ff) matches with the typical anisotropy 
structure found for round, non-rotating ellipticals by 
iGerhard et"aH fcOOll) . Along the m inor axis, NGC 4 807 is 
dominated by tangential motions. iGebhardt et alJ J2003T) 
report tangential anisotropy for some galaxies (mostly en- 
hanced azimuthal, but s u ppres s ed meridional dispersions) 
and a lso iDeionghe et al.| l|l996}). [Cretton fc van den Boschl 
J1998!) . IStatler. Deionghe fc Smecker-Hand Jl99Sl) 
ICappellari et alJ (120021). IVerolme et all (120021) 
and Icopm"" Cretton fc Emselleml H^Ool) find pre- 
dominantly tangential mo t ions. For NGC 3115 
lEmsellem. Deionghe fc Baconl lll999l) notice a z > ay. 
All these studies are based on data sets with different 
spatial sampling and use various dynamical modelling 
techniques, particularly differing in the amount and func- 
tional form of the applied regularization. As for the mass 
structure and the dark matter properties, a large and 
homogeneous sample of galaxies is needed to address the 
physical processes shaping elliptical galaxies. 
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